All-electron GW calculation based on the LAPW method: 
application to wurtzite ZnO 
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We present a new, all-electron implementation of the GW approximation and apply it to wurtzite 
ZnO. Eigenfunctions computed in the local-density approximation (LDA) by the full-potential lin- 
earized augmented-plane-wave (LAPW) or the linearized muffin-tin-orbital (LMTO) method sup- 
ply the input for generating the Green function G and the screened Coulomb interaction W. A 
mixed basis is used for the expansion of W, consisting of plane waves in the interstitial region and 
augmented-wavefunction products in the augmentation-sphere regions. The frequency-dependence 
of the dielectric function is computed within the random-phase approximation (RPA), without a 
plasmon-pole approximation. The Zn 3d orbitals are treated as valence states within the LDA; both 
core and valence states are included in the self-energy calculation. The calculated bandgap is smaller 
than experiment by ~1 eV, in contrast to previously reported GW results. Self-energy corrections 
are orbital-dependent, and push down the deep O 2s and Zn 3d levels by ~1 eV relative to the LDA. 
The d level shifts closer to experiment but the size of shift is underestimated, suggesting that the 
RPA overscreens localized states. 



I. INTRODUCTION 

Density-functional theory (DFT) provides a founda- 
tion for modern electronic-structure calculations, and the 
local-density approximation (LDA) is an efficient way to 
calculate the ground-state properties of material. How- 
ever, the LDA eigenvalues should not necessarily identi- 
fied with the quasiparticle (QP) energies, although eigen- 
value differences are often used to describe the excited 
state. Time-dependent DFT can in principle describe 
the excited state, but a good approximation for the time- 
dependent exchange-correlation kernel is not known. The 
GW approximation (GWA) of Hedin provides a prac- 
tical method to calculate the QP energy. Hybertsen and 
Louie presented the first GW calculation for real mate- 
rials in 1986 They employed eigenfunctions given by 
the LDA as input, using additionally a pseudopotential 
approximation. 

Several methods have since been developed within 
various band-structure-calculation schemes Calcu- 
lated QP energies typically agree well with experiment, 
for many kinds of materials. However, various kinds 
of approximations in addition to the GW approxima- 
tion itself are usually employed, whose adequacy has 
not been well tested. Here we present a method that 
makes minimal approximations in addition to the GW 
and random-phase approximations. We start with the 
LDA eigenfunctions generated by the full-potential lin- 
earized augmented-plane-wave (LAPW) 0] or a variant 
of the full-potential linearized muffin-tin-orbital (LMTO) 



method |j |6J , where in either case eigenfunctions are ex- 
panded in atomic-like local functions in the muffin-tin 
(MT) sphere regions, and in plane waves in the inter- 
stitial region. In order to treat the localized electrons 
accurately, the Coulomb interaction v and the screened 
Coulomb interaction W are expanded with a newly- 
developed mixed basis that consists of two kinds of ba- 
sis functions. One is the product-basis function devel- 
oped by Aryasetiawan and Gunnarsson , which is con- 
structed from the products of the local functions in the 
MT-sphere regions. The other is the interstitial plane 
wave (IPW) that takes zero in the MT-sphere regions 
and equals to the usual plane wave in the interstitial re- 
gion. With the mixed basis, we can treat localized elec- 
trons, even core electrons, on the same footing as the 
other extended electrons. The mixed basis is by con- 
struction essentially a complete basis for the expansion 
of W; therefore, for given eigenfunctions as input, our 
method can produce reasonably well converged QP ener- 
gies rather more efficiently than a method that expands 
W in plane waves alone. 

The LMTO method employs smooth Hankel functions 
as envelope functions ||, which are smoother and more 
accurate than the ordinary Hankel functions customary 
to the LMTO method. About 100 envelope functions 
were employed in this calculation. To check convergence, 
empty sphere was added to some of the calculations, 
mainly to enlarge the basis set. However, little differ- 
ence was found whether the empty sphere was added or 
not. 
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We apply this approach to wurtzite-type ZnO, whose 
valence bands consist of extended O 2p and Zn 4s or- 
bitals, and rather localized Zn 3d and O 2s orbitals. ZnO 
is important for optical device technology since the mate- 
rial is optically transparent and can be doped with both 
electrons and holes ||]. Compared with most other II- 
VI and III-V compounds such as ZnS, GaN, etc., the 
position of the cation d levels are rather high, and rela- 
tively close to the anion p-derived valence band maximum 
(VBM). The effect of the Zn 3d state is not negligible for 
the various properties. For example, the 3c? state couples 
to the VBM and pushes it upward, reducing the bandgap. 

This paper is organized as follows. Section II briefly 
describes the method. In Sec. Ill, we show LDA and 



GW results for ZnO, and compare them to experimental 
data and previously reported GW results. A summary is 



given in Sec. IV 



II. OVERVIEW OF THE GW CALCULATION 

The brief procedure of our GW calculation is presented 
here; it will be described in more detail elsewhere ||. In 
the Green's-function approach, the QP energy P kn and 
wavefunction f^ n (r) of a many-electron system are given 
as the solution of the equation 

[^kn-r-VH(r)]/ kn (r) 

- J S(r, r', P k „)/ k „(r')dV = 0, (1) 

where T is the kinetic-energy operator, Vh is the Hartree 
potential plus the electrostatic potential from nuclei, and 
£ is the self-energy. In the GW approximation, the self- 
energy is written as 

—e iu ' 5 G(r, r', u + w')W(r, r', u/). 

-oo 27T 

(2) 

It is convenient to divide £ into £ = £ x + £ c , with 
£ x = iGv the bare exchange term, and £ c = iGW c the 
correlation term. Here W c = W — v. 

We take a perturbative approach to find P kn . First, 
we solve the Kohn-Sham (KS) equation, 

[e kn —T — Vk(r) — K L c DA (r)]^ kn (r) = 0, (3) 

where V^ c (r) is the exchange-correlation potential in 
the LDA. We assume that £ — V^ DA is small enough 
to be treated a first-order quantity. Expanding £(P k „) 
around e kn , we obtain to first order 

Ek n ~ e k „ -I- ^ k „[(V' k n|E(e k „)|'!/' k r l ) ~ (V'knlV^f^lV'kn)], 

(4) 

where Z k „ is the renormalization factor defined by 

z kn = [i- (v>k„i(d£/Miu,=. k „iV'k„>r 1 - (5) 

We also estimate the first-order energy in the Hartree- 
Fock approximation (HFA) through 

£ k „ FA = £k„ + (V>k„|£ x |V>k„> - <^ k „|^ L c DA |^ k «>, (6) 



although the KS wavefunction may be different from 
the HFA one, and, therefore, E^ A is not the true self- 
consistent HFA value. 

In the augmented-waves methods, space is divided into 
the MT-sphere regions and the interstitial region. In 
both LAPW and the present LMTO method, the KS 
wavefunction is expanded as 

lm 13=1,2 G 

where the atomic- like local function A a i m p(r) is defined 
by 



baipi^Yimfr) in the a-atom MT-sphere 







otherwise 



(8) 

with orthogonal radial wavefunctions <t) a ip{r) (/3 = I or 
2), and spherical harmonics Y/ m (f). The IPW -Pq(i") is 
defined by 

{0 in the MT-sphere regions 

, . (9) 
e i( k +G) r j n ^ e interstitial region . 

The interactions, v and W, are well expressed by 
the product of two KS eigenfunctions in our per- 
turbative treatment. The product ^> kini (r)-0 k2 „ 2 (r) 
is expanded by the product of two local functions, 
4iim 1 /3 1 (r)A a i 2 m 2 ft(i'), in the MT-sphere regions, and 
by the product of two plane waves, Pq^ (r)P^ 2 2 (r), in the 
interstitial region. Following Aryasetiawan and Gunnars- 
son [ ?J the complete set of product functions is reduced 
by eliminating nearly linearly dependent ones. Taking 
the Bloch sum of the product functions, we finally ob- 
tain the product-basis function expressed as B^(r) in 
the a-atom MT-sphere region. The product of IPW's 
is also IPW in the interstitial region: (r)Pg 2 3 (r) = 

P Gi+G 2 ( r )- Thus, we obtain a mixed basis {AP^r)} = 
{B^(r), Pg(r)} which is suitable for expansions of v and 
W . The index i specifies a member of the basis and runs 
through G and afi. 

Because of the nonorthogonality of IPW's, the overlap 
integral of the mixed-basis functions, 

Sij = (M?\Mf) , 

is nonvanishing for i ^ j. We therefore define the dual- 
basis function: 



(10) 



The Coulomb interaction v is expanded as 

BZ 

k ij 

% (k) = (M?\v\Mf) . (II) 
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The Coulomb matrix w«(k) can be calculated by using 
the structure constants ||. 

The self-energy is calculated by using ek«, Vy(k) and 
(V'qnIV'q-kn' Mf) . The diagonal part of S x is given by 

(Vvl S *IVv) 

BZ occ. 

= - J] E H(V' q r l |V' q -kn'M^)^(k)(Af] i V q -k„'|V' q n)- 
k n' 

(12) 

The non-interacting polarization function Z) is expanded 
in the same manner as Eq.(|TT|): 

BZ occ. unocc. 

= EE E Wk„i<A k+qn' k+ q n' 

1 1 



k n 



V - Ck+qn' + e kn + ^ W + £k+qn' ~ £kn ~ ^ 

(13) 

We use the tetrahedron method for the Brillouin zone 
(BZ) summation in Eq.Q following Ref.|o|. The 



screened Coulomb interaction is given by W = (1 — 
vD)~ 1 v in the random-phase approximation (RPA). 
is also represented by the mixed basis. 

The correlation part of the self-energy is calculated as 

(Vvl^cMIVv) 

BZ All 
k n' ij 

Wg.(k,u/) 



2tt 



+ Ci/ - £q-kn' ± ^ 



dw', 



(14) 



where — i<5 is taken for the occupied states and +iS for the 
unoccupied states. We perform the frequency integration 
in Eq.(fl4]) with a method devised by Aryasetiawan ||. 

The evaluation of the exchange self-energy must be car- 
ried out carefully since the Coulomb interaction «y(k) 

shows a singular behavior u-y(k) oc — — as k — > 

0, where J7j(k) denotes the corresponding normalized 
eigenf unction. The singularity also exists in W. To 
handle the singular behavior, we use the offset T-point 
method || , which is essentially equivalent to the method 
to integrate the divergent part analytically ||. 

The QP energies are calculated with including the core 
contributions through the following equation, 



Ekn — £k)i + ^kn[ (V'knlS 



corcl 



^corc2+ valence 



\lpkn) 



(^kn|Kc DA (ntotal)|V^n> 



(15) 



where we divide the core states into two groups: corel 
is the deep core state, which affects the QP energies 
only through S x , and core2 is the relatively shallow core, 
which is treated on the same footing as the valence elec- 
trons. 



III. RESULTS FOR ZINC OXIDE 

We first describe the LDA calculation performed using 
the full-potential LAPW method. The local exchange- 
correlation functional of Vosko, Wilk and Nusair jll] is 
employed. The space group of wurtzite ZnO is PG^mc. 
The lattice constant and the MT-sphere radii are given 
in Table ffl. The angular-momentum in the spherical- wave 



expansion is truncated at Z n 



6 and L 



4 for 



the potential and the wavefunction, respectively. This 
l nm = 4 value is rather small, which gives errors com- 
pared with more accurate calculations with Z max = 7, but 
the differences of them are within 0.02 eV for the LDA 
band gap. The energy cut-off of the IPW is 16Ry for the 
wavefunction. We take 1152 fc-points in the first Bril- 
louin zone (BZ). The Zn (3d) 10 (4s) 2 and O (2s) 2 (2p) 4 
electrons are treated as valence electrons. 

The GW calculation is performed with 32 fc-points in 
the BZ. The energy cut-off of the IPW is 10 Ry for the 
Coulomb matrix. We treat 18 occupied bands and take 
into account 100 unoccupied bands. When producing the 
product basis, we ignore the product functions including 
<t>aif)=2 because the terms make small contributions. The 
Zn 3p states, which is relatively shallow in the core states, 
are chosen to be core2 electrons, i.e., the Zn 3p states are 
treated on the same footing as valence states and taken 
into account for the calculation of the correlation part 
of the self-energy. All the core and valence electrons are 
included into the calculation of the exchange part of the 
self-energy. In Sec. [II C, we check the convergence of the 



QP energies in fc-points, plane waves, unoccupied states, 
and product functions. 



A. LDA 

Unless otherwise stated, results in this section refer to 
the LAPW method. The LDA band structure for ZnO 
is shown in Fig. [l] and the density of states (DOS) is 
shown in Fig. |2j. All energies are given with respect to 
the top of the valence band. Around —17 eV, we obtain 
two bands originating from the O 2s states. The narrow 
bands between —6 and —4 eV consist mainly of the Zn 
3d orbitals, and the moderately dispersive bands from 
—4 to eV consist mainly of the O 2p orbitals. Fig. ^ 
shows significant p-d hybridization. The Zn 3d-derived 
bands are split into two groups, leading to a double-peak 
structure in the DOS. The lower peak is characterized by 
the strong p-d hybridization. The sharp upper peak be- 
tween —4.8 and —4.2 eV has strong Zn 3d character and 
the hybridization with the O 2p states is very small. The 
band gap opens between the 18th and 19th bands and 
the fundamental gap of 0.77 eV is located at the T point. 
The lowest two conduction bands consist mainly of the 
Zn 4s orbitals. The energy levels at the T point com- 
puted by the LAPW method are shown in Table |l|. Also 
shown are the levels computed by the LMTO method. 
Agreement is excellent; the bandgaps are nearly identi- 
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cal and the bands agree to within 0.1 eV over the entire 
valence bands. A further check was made using an en- 
tirely different kind of LMTO method E3|, and similar 
agreement was found. 

The band dispersion is similar to the result of the 
norm-conserving pseudopotential method by Schroer, 
Kriiger and Pollmann [|14| . One important discrepancy is 
the magnitude of the band gap: their fundamental gap is 
only 0.23 eV. Other reported LDA gaps are 0.93 eV cal- 
culated with the full-potential LAPW method Jl5| and 
1.15 eV within the LMTO- ASA method lUl. 



B. GWA 

The GW band structure computed with LAPW input 
is shown in Fig. || together with the LDA band structure. 
In the figure, all energies are given with respect to the top 
of the valence band (the valence-band top in the GWA 
was shifted up 0.49 eV relative to the LDA). We also 
show the LDA, GW and HFA energies at the T point in 
Table ||. The first and second bands are the O 2s bands, 
the 3- 12th are mostly Zn 3d character, bands 13-18 are 
of mostly O 2p character, and bands 19-20 the Zn 4s 
conduction bands. Because the Zn 3c? bands and the O 
2p bands overlap in energy, the characterization of the 
band is more complicated: e.g., at the F point, the LDA 
8th state has an O 2p character, while the 13th state has 
a Zn 3d character. 

The self-energy correction is sensitive to the character 
of the band, as seen from Table 0. The lowest O 2s 
bands are shifted by ~ — 1 eV; the Zn 3d bands by ~ — 1 
eV; the conduction Zn 4s bands by ~ +2 eV. The Zn 3d 
bandwidth shrinks, while the O 2p bandwidth and the Zn 
4s bandwidth are enhanced. The O 2s bandwidth does 
not change. 

The LAPW gaps are 0.77 eV, 2.44 eV and 11.39 eV in 
the LDA, GWA and HFA, respectively, as Table [n] shows, 
for the 32 fc-point mesh employed. The LMTO results are 
very similar. The experimental gap is 3.44 eV |fuj| . The 
large HFA gap reflects the neglect of screening the ex- 
change while the LDA gap is too small largely because 
it neglects the nonlocality in the (screened) exchange al- 
together @. The GW gap is closer to the experiment 
but still smaller by ~1 eV than the experimental gap. 
Previously reported GW gaps are somewhat larger. Us- 
ing a model GW approach, Massidda et al. found a gap 
of 4.23 eV |L5| and Oshikiri et al. found a gap of 4.28 
eV using an all-electron approach within the atomic 
spheres approximation. The major part of the error in 
our gap energy is thought to be derived from the overes- 
timate of the dielectric function, originating from under- 
estimated LDA gap used to generate it. We have found 
previously |5) that GW gaps are systematically under- 
estimated in semiconductors, with the error increasing 
with ionicity. Using an all-electron GW implementation 
within the PAW method (albeit including valence-only 
electrons in the GW calculation), Arnaud and Alouani 



p8| have also noted a tendency for the GW band gaps 
to be smaller than the experimental ones. For example, 
for Si, they find a gap of 1.00 eV. Using a fully converged 
512 fc-point mesh, our LAPW-GJU gap is 0.88 eV and 
our LMTO-GW gap is 0.89 eV. Wei Ku || has recently 
reported the LAPW-GIU gap of 0.85 eV using a 512 fc- 
point mesh. 



The centers of the O 2s band and the Zn 3c? band are 
-20.7 eV [Ha nd -8.81 eV |l|, respectively, in experi- 
ment. Table |I| shows that those GW bands are located 
higher by about 2 eV than the experimental positions. 
The HFA result is rather closer to the experiment. In 
our calculation, the dielectric function is overestimated 
for the reasons listed below. (1) In the RPA, the electric 
attractive force is neglected between excited electron and 
hole; there is no restoring force for the polarization, lead- 
ing to the overestimation of the polarization function. (2) 
The LDA band structure has a much smaller band gap 
than the experimental one. The electron-hole pair exci- 
tation energy is smaller than the real system. Since we 
use the LDA band structure to calculate the polariza- 
tion function, we overestimate it. (3) The Zn 3d and O 
2p bands are well separated in experiment, while they 
partly overlap in the LDA. The hybridization between 
the Zn 3d orbital and the O 2p orbital is overestimated. 
The hybridization between the Zn 4s orbital and the O 
2p orbital is also overestimated because of the smaller 
bandgap in the LDA. This property of the LDA wave- 
function makes the charge transfer between Zn and O, 
leading to the overestimation of the polarization func- 
tion. All those procedures of our calculation causes the 
overestimation of the dielectric function, and, therefore, 
the screened Coulomb interaction, W. In order to im- 
prove our calculation in the latter two points of (2) and 
(3), we need to perform a self-consistent calculation, that 
is, to use the QP energies and wavefunctions to calculate 
the screened Coulomb interaction W. 



In Table [II, we show the real part of £ at the 
T point. In the table, we take notations, S^ 1 ^ VA = 
(<MS(e k „)|Vk„> and S^ A = (^„|U x L c DA |Vk„>Lthe self- 

The 



energy corrections are given by Z^ n [£j^ WA 



normalization factor Z^ n is between 0.66 and 0.81, which 
is comparable to simple metals or semiconductors. As 
noted before, the self-energy corrections are negative for 
the valence bands and positive for the conduction bands 
and are larger for the localized states. Table |y| shows 
the decomposition of S^ VA at the V point to the core- 
exchange part ££™ rel , the exchange part ££„, and the 
correlation part The contributions of the core2 elec- 
trons are included in ££ n and The exchange part 
^kn nas a l ar g e discontinuity across the Fermi level, lead- 
ing to the wide gap in the HFA. The correlation part 
££ n are positive for the valence bands and negative for 
the conduction bands, leading to the reduction of the 
bandgap from the HFA value. 
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C. Convergence check 



The GW band structures described in Sec. Ill B have 



been calculated using 32 fc-points in the BZ, an IPW en- 
ergy cutoff of 10 Ry for the Coulomb matrix, 100 unoc- 
cupied states, and neglecting products including 4> a i}3=2- 
To check the convergence of the QP energies, we have 
performed GW calculations with some different condi- 
tions. The conditions we have used are 64 fc-points 144 
fc-points, the energy cut-off of 16 Ry for the IPW, or 200 
uno ccupied states. The results a re sh own in Table fv|- 
VII, compared to the result in Sec. Ill B. In any case, the 
errors are within 0.1 eV. We have also performed the cal- 
culation with the products including 4> a ip=2 for Zn. The 
improvement is quite small, within 0.005 eV. Our GW 
result thus shows good convergence; especially, a small 
number of IPW is enough to achieve a good result. 

We have c hecked the influence of the Zn 3p core states. 
In Sec.IIIB, we have treated the states as core2, i.e. they 
treated on the same footing as valence states. Now, we 
treat the Zn 3p states as corel (and include the products 
with 4> a if3=2 for Zn). The fundamental gap becomes 2.47 
eV and the Zn 3c? bands are shifted lower by ~ 0.1 eV. 
The effect of the Zn 3p electrons on the valence band 
through the correlation term in the self-energy is about 
0.1 eV. 



IV. SUMMARY 

We have presented a procedure for calculating the self- 
energy in the GWA with the mixed-basis expansion based 
on the full-potential LAPW and LMTO methods, and 



have applied the all-electron GW calculation to wurtzite 
ZnO. The mixed-basis method works well for this system 
which has both extended states and localized states; the 
GW calculation has a good convergence in various pa- 
rameters, and can be performed on a workstation-level 
computer. 

The GW bandgap of ZnO is smaller than experiment 
by ~ 1 eV. The self-energy correction is orbital dependent 
and the localized O 2s and Zn 3d states are lowered by ~ 
leV relative to the LDA values, while still higher than ex- 
periment. The GW calculation overestimates the screen- 
ing effect for localized states such as the Zn 3d states, 
because of the RPA and the LDA band structure. These 
errors are apparently systematic in zincblende semicon- 
ductors ||. We have indicated that a self-consistent cal- 
culation can improve the result. 

For a four-atom system having localized 3d orbitals 
such as wurtzite ZnO, we can complete the GW calcula- 
tion with 32 fc-points in the BZ within three days by 
using a DEC alpha 21264 667MHz workstation. Our 
GW code is available in a web site. The address is 



http://all.phys.sci.osaka-u.ac.jp 
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Es radius [a.u.] 




2.04 
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LDA GWA HFA 

Band n LAPW LMTO LAPW LMTO LAPW 



I 


—17.84 


— 17.74 


— 18.56 


—18.37 


—23.17 


2 


-17.06 


-16.97 


-17.78 


-17.67 


-22.42 




— <\ so 


— <\ 77 




— (S 

\J.O-L 




5 


—5.74 


—5.70 


—6.58 


—6.49 


—9.99 


( i 7 
u. / 


—5.62 


—5.59 


—6.46 


—6.34 


—9.80 


8 


-5.53 


-5.61 


-5.65 


-5.59 


-5.58 


9,10 


-4.70 


-4.63 


-5.90 


-5.82 


-9.87 


11,12 


—4.48 


—4.42 


c; co 

tJ.UO 


— ^ fil 


— Q 7<\ 


13 


-4.19 


-4.12 


-5.20 


-5.13 


-8.92 


14,15 


-0.74 


-0.79 


-0.84 


-0.87 


-0.89 


16 


-0.10 


-0.09 


-0.08 


-0.03 


-0.00 


17,18 

















19 


0.77 


0.78 


2.44 


2.44 


11.39 


20 


5.13 


5.12 


7.19 


7.25 


17.13 



TABLE II: The LDA, GW and HFA energies (in eV) at k = 
(0, 0, 0). All energies are given with respect to the top of the 
valence band. The LAPW-GW and LMTO-GW results use 
32 a nd 2 16 fc-points, respectively. As shown in Tables [v|, |vi| 
and [Vll|, the LAPW-GW fundamental gap would be 0.05 to 
0.1 eV smaller if the combined effect of complete convergence 
in fc-points, IPW and number of unoccupied states were taken 
into account. 
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TABLE III: GW self-energies S£^ A at k = (0,0,0), to- 
gether with the LDA exchange-correlation term £^° A and 
the renormalization factor Z^ n . The corrections are given by 

rv r v GWA v LDAi 
^knl^kn — ^kn J- 



Band n 


S« n WA [eV] 


[eV] 




Corrections 


i 
i 


99 HQ 


90 


u.oy 


1 99 
— l.ZZ 


o 
L 


93 98 


91 zi^ 


U.DD 


1 91 


3,4 


-35.25 


-33.37 


0.70 


-1.32 


5 


-35.74 


-33.83 


0.70 


-1.33 


D, < 


Od. ( 4: 


33 RA 

00.04: 


U. lU 


— l.OO 


Q 

o 


1 Q «0 
— lo.DZ 


1 1 1 A 


U. f U 


— U.ul 


9,10 


-41.70 


-39.36 


0.72 


-1.69 


11,12 


-42.46 


-40.11 


0.72 


-1.69 


13 


-40.66 


-38.55 


0.71 


-1.50 


14,15 


-27.44 


-26.65 


0.75 


-0.59 


16 


-27.93 


-27.30 


0.75 


-0.47 


17,18 


-28.79 


-28.13 


0.75 


-0.49 


19 


-12.87 


-14.32 


0.81 


1.18 


20 


-12.45 


-14.39 


0.81 


1.57 



TABLE IV: Core-exchange part E£™ rel , exchange part E£ n 
and correlation part E£. n of the GW self-energy at k = 
(0, 0, 0). The contributions of the core2 electrons are included 
in E£ n and E c kn . 



Band n 


££T cl [eV] 


SL [eV] 


£L [eV] 


1 


-1.85 


-29.84 


9.60 


2 


-2.02 


-30.81 


9.55 


3,4 


-5.04 


-38.51 


8.29 


5 


-5.17 


-38.93 


8.36 


6,7 


-5.18 


-38.87 


8.31 


8 


-1.24 


-22.57 


5.20 


9,10 


-6.83 


-43.73 


8.85 


11,12 


-7.01 


-44.39 


8.95 


13 


-6.55 


-42.75 


8.65 


14,15 


-2.85 


-29.97 


5.39 


16 


-2.92 


-30.31 


5.30 


17,18 


-3.14 


-31.02 


5.37 


19 


-1.45 


-8.27 


-3.15 


20 


-1.44 


-6.97 


-4.04 
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TABLE V: Dependence on the QP energies on the number 
of fc-points, for representative states at the P point. The cal- 
culation with 32 k points is the same one in Table & The 
LMTO method was used to compute QP energies for finer 
meshes, up to 216 fc-points (6x6x6 mesh). The LMTO-GW 
gap changed by —0.1 eV going from 64 to 216 fc-points; as- 
suming the same convergence in fc-points within the LAPW 
method, we estimate the converged LAPW-G14 7 fundamental 
gap to be 2.32 eV, keeping all other parameters fixed. The 
LAPW-GW and LMTO-GW energies agree with each other 
to within 0.1 eV, with the error approximately tracking dif- 
ferences in the LDA eigenvalues. 



Band n 144 k 64k 32 k 

1 -18.51 -18.54 -18.56 

3 -6.60 -6.60 -6.62 

8 -5.64 -5.62 -5.65 

14 -0.84 -0.83 -0.84 

18 

19 2.35 2.42 2.44 



TABLE VI: Dependence of the QP energies on the number of 
IPW used to make the the Coulomb matrix, for representative 
states at the P point. Data with the 10 Ry cutoff is the same 
as in Table |j| 



Band n 16 Ry 10 Ry 

1 -18.55 -18.56 

3 -6.59 -6.62 

8 -5.64 -5.65 

14 -0.84 -0.84 

18 

19 2.47 2.44 



TABLE VII: Dependence of the QP energies on the number of 
the unoccupied states, for representative states at the P point. 
The calculation with 100 unoccupied states is the same as in 
Table p. 



Band n 200 unocc. 100 unocc. 

1 -18.56 -18.56 

3 -6.54 -6.62 

8 -5.62 -5.65 

14 -0.83 -0.84 

18 

19 2.45 2.44 



LDA 
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FIG. 1: LDA band structure in wurtzite ZnO. 
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FIG. 2: Density of states in wurtzite ZnO. 




FIG. 3: The GW band structure (solid lines) and the LDA 
band structure (dashed lines) for wurtzite ZnO. The self- 
energy corrections are calculated at the points marked '*'. 



